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ABSTRACT 



This paper evaluates five maneuvering strategies for an 
evader in a two dimensional continuous evasive game. An 
evader's movement, in two dimensions at constant speed, is 
simulated by choosing courses that are independent and dis- 
tributed according to one of the five course change rules . 
The times between course changes are independent exponential 
random variables. An improvement over a previously estab- 
lished least upper bound on evader survivability is pre- 
sented. The bulk of this paper is the presentation of 
evader survivability as a result of the maneuvering strategy 
he employs. 
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I. 



INTRODUCTION 



A classic military problem with current application is: 
how to maneuver a mobile target in order to prevent success- 
ful prediction of its position. For example how should an 
aircraft carrier maneuver to maximize its survivability 
against an ICBM attack? 

Of primary importance to this problem is the time lag 
between the decision to fire at a target, made at time t, 
and the warhead detonation at time t+TL. This lag is called 
time-late and denoted as TL. The decision to initiate such 
an action is based in part upon knowing the target's posi- 
tion at time t while the effectiveness of that action is 
dependent upon the error in predicting that position at t+TL. 

A. ICBM VERSUS AIRCRAFT CARRIER 

Consider the ICMB versus aircraft carrier problem. As- 
sume that the carrier operates in mid-ocean, possesses great 
endurance, can make sharp turns and has as its only kinematic 
restriction a constant speed, v. The carrier is referred to 
as the evader, E. The attacker, P, continuously observes E 
from a nearby unarmed trawler. P's weapon is a land based 
ICBM which has perfect accuracy, produces a lethal area A 
(this lethal area can have any shape) but does not have a 
mid-course guidance nor homing capability. Although P knows 
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E's position past and present, he also knows that his ICBM 
will have a one hour time-late.^ 

P's problem is to predict E's position at time t+TL in 
order to initiate an attack at time t. E's problem, since 
he will not know of an attack until he observes the deto- 
nation, is to maneuver in such a way to confound P's predic- 
tion. E's movement should be random because any nonrandom 
movement would be vulnerable to extrapolation. 

B. AN UPPER BOUND ON EVADER SURVIVABILITY 

A lower bound on P's success, p^,, and consequently an 
upper bound on E's survivability, 1-p^/ is known. P knows 
E's position at t+TL will be interior to or on a circle of 
radius v*TL centered on E's position at t (which is known). 
The enemy can then gaurantee a kill probability of at least 
A/ir(v*TL) by choosing, a lethal area randomly within the 
circle of uncertainty (a wedge of random orientation will 
do) [Ref. 4] . Then the maximum survivability E could attain 
would be 1 -A/tt (v*TL) ^ . 

This upper bound on E's survivability is a function only 

of P's strength, that is the magnitude of A (the value 
2 

ir(v*TL) is constant). Figure 1 is a graph of that survival 
function versus P's probability of success. 

■^Time-late is a summand of many factors in the command and 
control problem as well as the missle flight time. For the 
purpose used here it is sufficient to treat time-late in 
total and not consider its decomposition. 

2 

If E knows the time at which an attack is initiated then 
the problem is trivial. 
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Survivability Versus Attacker Strength 

FIGURE 1 

Unfortunately the specific set of rules for 
by to achieve the survivability depicted in Fig. 
known . 




to maneuver 
1 is not 
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II. BACKGROUND 



A. DISCRETE EVASIVE GAMES 

It was felt that a two dimensional continuous game was 
too difficult to solve. So the approach was to solve a 
similar but easier game. It was assumed that the ocean was 
a linear set of discrete points and that E's mobility con- 
sisted of being able to jump either right or left to an 
adjacent point [Refs. 2 and 3]. Time-late became an integer 
number of jumps that E could make between the attack de- 
cision and warhead detonation. As in the exact problem E 
had no knowledge of the decision to attack. 

Discrete problems of this type are called evasive games 
and are classified by the number of jumps constituting time- 
late. A one step discrete evasive game means E can move 
either left or right one jump prior to detonation. Game 
theory provides an immediate solution to that game; the 
value is pj,=0.5. Also there is an optimal strategy for E 
in that game. At each jump E should go left with probability 
one-half or right with probability one-half. Employment of 
such a rule for each jump guarantees E a survival probability 
of at least the game value, regardless of P's firing rule. 

Having easily solved the one step game interest was 
focused on the two step game. In this game E was allowed 
two jumps during time-late. The solution has been obtained, 
but it was not as easy to achieve as that of the one step 
game [Ref . 2] . 
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The direction of the analysis was clear. Knowledge 
gained from solving the more simple games would be a step- 
ping stone to the solution of the more difficult games . 
Eventually the assumptions of linearity and discreteness 
could be relaxed. Progress, however, in solving the dis- 
crete games has been very slow. Researchers are presently 
embroiled in solving the three step game. 

B. A CONTINUOUS EVASIVE GAME 

Washburn in Ref. 1 presented a different approach to 
analyzing the ICBM versus aircraft carrier problem. He de- 
veloped the probability density function of a particle 
moving continuously in a two dimensional medium subject to 
a specific set of maneuvering rules . That approach was 
distinct from previous work because it addressed the exact 
problem. The results presented in this paper are an ex- 
tension of that approach. 

Washburn proposed the following strategy for E and 
analyzed the evader's subsequent survivability, first as a 
function of the attacker's strength and secondly as a 
function of a specific strategy parameter. E was to select 
courses from a uniform distribution. Course changes were 
to occur as a Poisson process. This meant the time between 
course changes would be exponentially distributed. 

The exponential time between course changes has an 
intuitive appeal because of the memoryless property of that 
distribution. The memory less property is: the probability 
that E will not change course in the future, given that he 
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has not changed course for some observation period, is in- 
dependent of the length of that period. Therefore, use of 
the exponential time distribution should serve to confound a 
prediction of future position by extrapolation, regardless 
of the course distribution used. For example, suppose P's 
attack decision rule required that he observe E maintain a 
constant course for at least five hours prior to initiating 
an attack. Such a procedure would not improve P's prob- 
ability of a kill because the only information of any benefit 
at the time an attack is initiated is E's position and last 
course. 

Figure 2 is a graph of the evader's maximum attainable 

survivability when using the strategy of uniform courses 

and exponential times versus P's strength for the optimal 

3 

exponential parameter, A. For comparison the upper bound 
on survivability is also presented. 

"If we assume the parameters A,v,TL are known to both 
sides, then the evader can select A to maximize the sur- 
vivability. The evader will clearly be in trouble if he 
makes A too small, because the kill probability is at least 
exp(-A*TL).^ On the other hand making course changes too 
frequently will lead to a density function that is highly 
peaked at the origin, which is equally undersirable . . . " 

[Ref. 1]. 

3 

The exponential parameter, A, is the inverse of the mean 
time between course changes . 

4" 

Exp(-A*TL) is the probability that E will make no course 
change during TL. 
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X -Pk 




Minimum p k for Optimal Parameter A Versus Attacker 
Strength; Strategy Is Uniform Courses 
and exp (A) Times 

FIGURE 2 



"It is not known whether or not the optimal strategy for 
E is a Poisson strategy of the type just considered, or even 
whether the uniform distribution on angles is optimal within 
the class." [Ref. 2] 
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III. RESULTS OF SIMULATING THE CONTINUOUS EVASIVE GAME 



The goal of this research was to determine if the uni- 
form course rule was optimal within the class of Poisson 
strategies or to find a better rule if it was not. For 
another rule to be better than the uniform one it would 
have to enable E to attain a higher survivability than that 
shown in Fig. 2. 

The methodology of solving discrete games had not ad- 
vanced sufficiently to achieve those goals. The methods 
used by Washburn were also not suitable because of the 
mathematical difficulty of the problem. 

The Poisson class of strategies is, however, uniquely 
suited for analysis by computer simulation. This is due to 
the memory less property of the exponential times between 
course changes. The results of such a simulation are pre- 
sented in this paper in the form of improved survivability 
for the evader. 

Five course change rules were evaluated in the simula- 
tion. The probability density function of each is pre- 
sented on subsequent pages along with a graph of the re- 
sulting attacker's p^ versus attacker strength, denoted as 
S=A/tt (v*TL) ^ These graphs are for representative values 
of the parameter a. Alpha, a, is the product of time-late 
and the exponential parameter A. 

5 

Course changes were assumed to be independent of the 
underlying Poisson process. 
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The simulation established a new lower bound on the 
the evader could yield to the attacker. Figure 3 is a graph 
illustrating the improvement of the new bound over that 
found by Washburn. This new bound is the least lower bound 
of all the p k versus S curves for the four values of a and 
all rules simulated. Examples of such curves are shown in 
Figs. 5-9. The old bound in Fig. 3 is the lower bound of 
the p^ curves simulated using the uniform rule exclusively. 

The minimal p^ as a function of a, over the five rules, 
is shown in Fig. 4 for three specific attacker strength 
levels. The curves in Fig. 4 confirm the statement in Ref. 

2 that the evader should turn most often, that is a should 
be highest, against the weakest opponent. 

Comparisons of the pj. curves of the different rules 
have shown that the uniform rule is not optimal for all 
attacker strengths and a values. For example, the reverse 
course rule (see Fig. 7) was shown to be better than the 
uniform rule for ot=l. The difference between these two 
rules is illustrated in Fig. 10. Against a weak opponent 
(strength ^.6) the evader would do better to use the re- 
verse course rule rather than the uniform rule for a=l. 

This improved survivability occurs because the reverse 
course rule created a "flatter" evader position density. 

None of the five course rules evaluated produced a con- 
sistently smaller p^ for all values of attacker strengths 
and a. However, the left-right rule was dominated by the 
o£her four in all cases. Each of the five course rules was 
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simulated at six different levels of a by varying ET . The 
graphs of p^. versus A/ir(v*TL)^ for each simulation are pre- 
sented in the Computer Output section along with tabulated 
p^'s for specific attacker strength levels as a function 
of the course rule and the a used. 




New Bound on p^. That the Evader Can Yield 
as a Function of Attacker Strength 

FIGURE 3 
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1.0 n 




FIGURE 4 
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Density Function of Uniform Rule 
FIGURE 5A 




FIGURE 5B 




Density Function of Modified Left-Right Rule 



FIGURE 6 A 



Pk 




FIGURE 6B 
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Density Function of Reverse Course Rule 



FIGURE 7 A 



Pk 




FIGURE 7B 
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Probability [Course Change = tt/2] = 0.5 
Probability [Course Change = 3tt/2] = 0.5 



Density Function of Left-Right Rule 



FIGURE 8A 



Pk 




FIGURE 8B 
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Density Function of Truncated Uniform Rule 



FIGURE 9A 




FIGURE 9B 
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Pk 




A Comparison of the Tradeoffs Between the Uniform 
and Reverse Course Rules As a Function of Attacker Strength 



FIGURE 10 
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IV. SIMULATION MODEL 



Let E's position at time t be at X=0 , Y=0 of a cartesian 
coordinate system. Then rotate the axis until the positive 
abscissa is aligned on and pointing in E's last direction 
of travel. Then regardless of the course rule used by E 
there is at least exp (-a) probability that he will be at 
X=v*TL, Y=0 at time t+TL. This is the probability that E 
does not change course during TL. So the attacker's prob- 
ability of killing E by centering a small portion of his 
lethal area at (v*TL,0) is exp (-a). Then add to exp (-a) 
the integral of E's density function over that part of the 
uncertainty circle where the rest of the lethal area is 
targeted to determine the total p^,. The total amount of 
lethal area allocated in this manner is the numerator, A, 
of the attacker's strength function S. 

A computer program to simulate the above computation 
was written in FORTRAN and run on the IBM 360. The program 
consisted of two phases, first the play of a strategy and 
secondly the scoring of that play. 

A. SIMULATION OF A STRATEGY 

To create a single play the evader's track was simulated 
from an initial position at time t to the resultant position 
at t+TL. The track was the result of a specific maneuvering 
strategy being simulated. A strategy was made up of two 
decision rules. The first rule determined the times be- 
tween course changes which were exponential random variables 
with mean ET=1/A. This rule was common throughout all of 
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the simulations, although the parameter ET was a variable. 
The second rule determined the magnitude of subsequent 
course changes. The only other kinematic restriction was 
a constant speed, v, for the evader. 

To start a play datum was initialized by setting t=0, 
X(0)=0, Y(0)=0 and an initial course, c^ , was selected from 
the uniform distribution. Then a sequence of exponential 

g 

times [t^]^_^ n(TL) were generated and a sequence of course 
changes, [ 3 ^_2 n(tL)' were 9 enerat ®d using the course 
change rule. The evader's position at t=TL was 
N-l 

X (TL) = Y v* t^ • (Cos (c^) -Cos (c N ) ) + v*TL*Cos(c n ) 1.1 
i=l 



N-l 



Y(TL) = Y v* t. • (Sin (c • ) -Sin (c N ) ) + v*TL*Sin(c N ) 1.2 
i=l 1 



B. SCORING THE PLAY 

Once E's position at t=TL was determined that observa- 
tion was scored. To score a play the coordinates of E's 
position, X (TL) , Y (TL) , had to be transformed. The purpose 
of the transformation was to make the observed position in- 
dependent of the particular initial cource c-^. The trans- 
formation was a rotation of the coordinate axis about the 
datum so the positive abscissa would be aligned in the 
direction of the initial course. The transformed position 
was : 

^The number of elements of this sequence, N(TL) , is a 
Poisson random variable. 
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2.1 



X' (TL) = X(TL)Cos(c ]L ) + Y(TL)Sin(c ) 

Y' (TL) = Y (TL)Cos (c 1 ) - X(TL)Sin(c ) 2.2 

An example of this transformation is shown in Fig. 11. 

If E had not made a course change during TL then his 
transformed position would have been 

X' (TL) = v* TL 3.1 

Y' (TL) =0 3.2 

. 7 

A grid system of square cells was placed over the 
playing area of Fig. 11B and a determination was made as to 




ORIGIN: Datum 

B: X (TL) , Y (TL) 

C^: Initial Course 

SAMPLE PLAY PRIOR 
TO TRANSFORMATION 




ORIGIN: Datum 

B: X' (TL) , Y' (TL) 

A: (v*TL,0) 

SAMPLE PLAY AFTER 
THE TRANSFORMATION 



FIGURE 11A FIGURE 11B 

which cell X' (TL) , Y' (TL) was in. Each cell of the grid 
had an associated value which would represent the number of 
times a play resulted in an observation in that cell. When 

7 

The cumulative cell area was normalized by the factor 
it (v*TL) 2 . 
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the appropriate cell was determined, for the play being 
scored, that value was incremented by one. 

A simulation run was composed of 16810 plays and scoring 
iterations for a strategy utilizing a specific course change 
rule. The input variables for a run, besides the course 
rule were v, TL and ET . At the end of a simulation run the 
grid system was a two dimensional histogram of E's position. 
The frequencies in the histogram were then ordered, accumu- 
lated and normalized to achieve cumulative cell probabilities. 
The ordering corresponded to the conservative assumption that 
P could divide his lethal area and target only those cells 
with the higher probabilities. These probabilities were 
then plotted against the cumlative cell area they repre- 

o 

sented. This graph was labeled pj, versus attacker strength 
and was the primary output of the program. 

Thirty simulation runs were made to investigate five 
different course change rules and six different values of 
ET. For all runs the following constant values were main- 
tained, v=5 knots and TL=2 hours. The graphs from those 
runs are included in the Computer Output section. Each 
course change rule was simulated many different times, the 
only difference between runs being the pseudo random num- 
ber generator seeds. This was done to check for the vari- 
ability of the p^ graph for that rule. In all such runs 



The amount of area in a cell was denoted as cell size and 
equal to (2*TL/41)2. tl and v were held constant there- 
fore cell size was always approximately 1/4 nm^. 
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the resultant graphs were so similar that any difference 
was not distinguishable. For simulation run required two 
minutes and forty-five seconds using 125 K on the IBM 360. 
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V. CONCLUSION 



It was determined that the uniform course change rule 
is not optimal within the Poisson class of strategies . Com- 
parison of various rules, see Fig. 10, showed that under 
certain conditions a single rule such as the reverse course 
change rule is better. Also the new lower bound on pj, , 
taken over all rules and values of a simulated, is an im- 
provement on the lower bound achieved from only the uniform 
rule. 

Certainly not all the possible course change rules were 
simulated. The rules evaluated, however, were representa- 
tive of the broad class of possible rules. When compared 
to the uniform rule, all but one of the other rules showed 
that E could improve his situation if he knew P's strength 
by selecting the better rule for that encounter. The one 
rule that was consistently dominated was the left-right 
rule . 

It remains unknown whether or not the optimal strategy 
is a Poisson strategy. 
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APPENDIX A 



FLOW DIAGRAM 



START 



IX -- 55555| 



■pseudo-random number seed, an input variable 

t 1 1 — 



| V — 5 . 0 \ | evader speed in knots \ 



| TL -- 2.0 1 time- late in hours | 



{ ET -- l7o | 






PP (41,41) 
PP(i,j) — 0 
for all i,j 



•Jmean time between course changes, an input 
jvariable 



Jplaying area and histogram; each cell is 

•a square with dimension (2*V*TL)/4l 



INDEX(41,41) 
INDEX(i , j ) -- 0 
for all i,i 


1 INDEX is a book-keeping record of PP cell 

Ifeasibility 
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0 - 



l c -- 


7| 


C + C 1 — 


4 




Generate 


a course 


change 


c by the 


rule being sim- 


lated 





Generate cv~U(0,2 ) 


} 


t 


Theta — c I 




t 


> 


Generate t^Exp(l/ET) 


> 


L 


T — min(TR, t) 


_ > 


f 


Clock — 


■ Clock + t 



x — x + (V*TL)*Cos(c) 
y — y + (V*TL)*Sin(c) 



|TR — TL-Clock 


* No 


Is Clocl 


< = TL 



.Yes 



x — x*Cos (Theta)+y*Sin (Theta) 
y — y*Cos (Theta)-x*Sin(Theta) 

1 



Determine which column of PP the x co- 
ordinate is in, e.g. column 1 



Determine which row of PP the y coordi- 
nate is in, e.g. row m 



PP (m, 1) — PP (m, 1) + 1 



z 



Q. 



Sort and order PP 



r 



Sum PP 



STOP «- 



x t 

\ Graph: versus attacker strength, \ 

\ i.e. p^ versus normalized lethal area \ 
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COMPUTER OUTPUT 
TL/ET 



Rule 


1.0 


2.0 


3.0 


4.0 


Uniform 


.729 


.580 


.501 


.486 


Modified 

Left-Right 


.737 


.580 


.501 


.486 


Reverse 

Course 


.690 


.572 


.548 


.564 


Lef t-Right 


.917 


.768 


.650 


.588 


Truncated 

Uniform 


.697 


.556 


.509 


.517 


column 

minimum 


.690 .556 .501 

for attacker strength equal 0.2 
Table I 
TL/ET 


.486 


Rule 


1.0 


2.0 


3.0 


4.0 


Uniform 


.854 


.784 


.752 


.760 


Modified 

Left-Right 


.878 


.784 


.752 


.768 


Reverse 

Course 


.838 


.791 


.799 


.831 


Left-Right 


.964 


.909 


.870 


.854 


Truncated 

Uniform 


.846 


.768 


.768 


.815 


column 

minimum 


.838 .768 .752 

for attacker strength equal 0.4 


.760 



Table II 
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TL/ET 



Rule 


1.0 


2.0 


3.0 


4.0 


Uniform 


.933 


.909 


.901 


.917 



Modified 

Left-Right 


.948 


.909 


.901 


.917 


Reverse 


.933 


.917 


.925 


.948 


Course 










Left-Right 


.987 


.987 


.980 


.987 


Truncated 


.933 


.909 


.909 


.933 


Uniform 










column 


.933 


.909 


.901 


.917 



minimum 



Pk for attacker strength equal 0.6 
Table III 



36 



Pk 




Pk versus Attacker Strength 
Course Change Rule: Uniform 



ET = 2.0 



37 



Pk 




A/7T (v -TL) 2 



pj. versus Attacker Strength 
Course Change Rule: Uniform 



ET = 1.33 
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Pk 

1.0 




Pk versus Attacker Strength 
Course Change Rule: Uniform 



ET = 1.0 



39 



?k 




A/tt (v -TL) 2 



Pk versus Attacker Strength 
Course Change Rule: Uniform 



ET = 0.8 
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Pk 




A/tt (v *tl) 2 



versus Attacker Strength 

Course Change Rule: Uniform ET = 0.66 
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Pk 




Pk versus Attacker Strength 

Course Change Rule: Uniform ET = 0.5 
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A/ Tt (v-TL) 2 



Pk versus Attacker Strength 
Course Change Rule: Modified Left-Right 



ET = 2.0 



43 



Pk 




Pk versus Attacker Strength 

Course Change Rule: Modified Left-Right ET = 1.33 
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Pk 




Pk versus Attacker Strength 

Course Change Rule: Modified Left-Right ET = 1.0 
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i : 





versus Attacker Strength 

Course Change Rule: Modified Left-Right ET = 0.8 
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A/tt (v-TL) 2 

versus Attacker Strength 

Course Change Rule: Modified Left-Right ET = 0.66 
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versus Attacker Strength 
Course Change Rule: Modified Left-Right 



ET = 0.5 
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Pk versus Attacker Strength 
Course Change Rule: Reverse Course 
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2 

A/ 71 (v*TL) 

versus Attacker Strength 

Course Change Rule: Reverse Course ET = 1.33 
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Pk 

1.0 - 




0.8 - 


/ / 


0.6 - 


/ / 


0.4 - 


/ / 


0.2 - 


/ / 
t / 

V . l _ 1 1 i i 




0.2 0.4 0.6 0.8 1.0 

A/Tv (v • TL) 2 

Pk versus Attacker Strength 



Course Change Rule: Reverse Course ET = 1.0 
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Pk 

1.0 



0.8 



0.6 



0.4 



0.2 




0.2 0.4 0.6 0.8 1.0 



A/TC (v-TL) 2 



versus Attacker Strength 
Course Change Rule: Reverse Course 



ET = 0.8 
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A/7C (v-TL) 2 



versus Attacker Strength 
Course Change Rule: Reverse Course 



ET = 0.66 



53 



Pk 

1.0 



0.8 * 



0.6 - 



0.4 - 



0.2 




0.2 



0.4 



0.6 



0.8 



1.0 



A/ it (v-TL) 2 



versus Attacker Strength 
Course Change Rule: Reverse Course 



ET = 0.5 
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A/ 75 (v-TL) 2 

p k versus Attacker Strength 

Course Change Rule: Left-Right ET = 2.0 
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A/ Tt (v*TL) 2 



PIq versus Attacker Strength 
Course Change Rule: Left-Right 



ET = 1.33 



56 



Pk 




A/ 71 (v*TL) 2 



versus Attacker Strength 
Course Change Rule: Left-Right 



ET = 1.0 



57 




versus Attacker Strength 

Course Change Rule: Left-Right ET = 0.8 
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A/ It (v*TL) 2 . 



versus Attacker Strength 
Course Change Rule: Left-Right 



ET = 0.66 
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Pk 

1.0 




p, versus Attacker Strength 
k 

Course Change Rule: Left-Right 



ET - 0.5 
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A/n (v-TL) 2 

versus Attacker Strength 

Course Change Rule: Truncated Uniform ET = 2.0 
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versus Attacker Strength 

Course Change Rule: Truncated Uniform ET = 1,33 
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A/ 7\ (v-TL) 2 

versus Attacker Strength 

Course Change Rule: Truncated Uniform ET = 1.0 
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A/lx (v-TL) 2 

Pk versus Attacker Strength 

Course Change Rule: Truncated Uniform ET = 0.8 
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Pk 




versus Attacker Strength 
Course Change Rule: Truncated Uniform 



ET = 0.66 
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A/7T (v-TL) 2 

versus Attacker Strength 

Course Change Rule: Truncated Uniform ET = 0.5 
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COMPUTER PROGRAM 



A SIMULATION OF THE CONTINUOUS EVASIVE 



GAME FOR THE POISSON CLASS OF STRATEGIES 



INTEGER * 2 PP, INDEX 
INTEGER * 2 CUMSUM 
INTEGER * 2 NCOUNT 
INTEGER *2 XP,QQ 
INTEGER * 4 CRITE3 
INTEGER * 4 TOTCEL 
REAL * 4 IN 
REAL LABEL/’ ’/ 

R£AL*8 I TI TL p ( 12)/'CL0THI£R' » ' THOMAS'.' J, ',3*' 
X ',' STRATEGY' , ' TRNCRS • , ' ( LAMDA ' , ' * TL )= ', 

X2*' 4.0'/ 

BY MAKING THE HISTOGRAM INTEGER STAR TWO THE AMOUNT 
OF STORAGE REQUIRED FOR THE PROGRAM WILL BE REDUCED. 



DIMENSION PP (4 1 ,41 ), INDEX (41,41 ) 

DIMENSION XP(1681) ,00(41 ,41) 

EQUIVALENCE ( P P ( 1 , 1 ) , XP ( 1 ) ) 

DIMENSION EEX ( 2 ) , EEY ( 2 ) , EX ( 692 ) , EY ( 692 ) , EX X< 692 ) 
DIMENSION E YY( 692) 



PP IS THE PLAYING AREA AND HISTOGRAM. EACH CELL IN 
PP REPRESENTS AN AREA OF SIZE CELSIZ. INDEX IS A 
BOOK-KEEPING ARRAY USED TO DETERMINE WHICH CELLS IN 
PP ARE FEASIBLE, E.G. IF INDEX(I,J)=0 THEN PP(I.J) IS 
AN INFEASIBLE CELL. SIMILARLY IF INDEX(I,J)=1 OR 2 
THEN PP(I,J) IS A FEASIBLE CELL. 

PP IS A SQUARE WITH THE TIME-LATE CIRCLE INSCRIBED IN 
IT. THEREFORE PP HAS SOME CELLS THAT ARE NOT IN SIDE 
THE TIME-LATE CIRCLE AND FOR THOSE CELLS THEIR INDEX 
VALUE IS ZERO. FOR THOSE CELLS IN PP THAT ARE CON- 
TAINED ENTIRELY- IN THE CIRCLE THE ASSOCIATED INDEX 
VALUE IS ONE. AND LASTLY THERE ARE THE CELLS IN RP 
THAT ARE ON THE EDGE OF THE CIRCLE, THESE ARE CALLED 
EDGE CELLS AND HAVE THE INDEX VALUE OF TWO. 



NC CUNT =0 
PI = 3. 1415926 

DATA IX,V,TL, ET/37915,05 .0,02.00,00.5/ 

V IS PARTI CLe SPEED 
TL IS THE TIME-LATE 

5T IS THE MEAN TIME BETWEEN COURSE CHANGES 

IX IS THE SEED FOR THE PSEUDO-RANDOM NUMBER GENERATORS 



DFLTAX=(2.C*V*TL)/41.0 

DELTAX IS THE HORIZONTAL DIMENSION OF A CELL. 
C=LSIZ=DELTAX**2 

CELSIZ IS THE AMOUNT OF REAL AREA CONTAINED IN A CELL. 
K = 20 

THE NUMBER OF CELLS IN EITHER DIMENSION IS ALWAYS A 
CONSTANT 41. THIS CONSTANT IS EXPRESSED AS AN ODD 
INTEGER AND WRITTEN IN THE FOLLOWING MANNER, 

41 = 2 * ( K ) + 1. THEREFORE K = 20. 



ZERO OUT THE HISTOGRAM 
DO 5 1=1,41 

DO 4 J=1 , 41 

INDEX ( I, J ) = 0 
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PP(I, J)=0 
00(1 ,J) =0 

4 CONTINUE 

5 CONTINUE 



NOW TO CHECK EACH OF THE CELLS IN THE HISTOGRAM AND 
DETERMINE ITS ASSOCIATED INDEX VALUE. 



QUADRANT I 
DO 42 J = 2 1 , 41 
DO 41 1=1,21 

0UT=( V*TL-J ;i DELTAX)**2+( ( I -1 ) * DE LT AX -V*T L ) **2 
IN =(V*TL-(J-1 )*DELTAX)**2+( I*DELTAX-V*TL ) **2 
RADSQ=( V*TL)**2 

IF( ( OUT. LE. RAD SO). AND. (IN. LE.RADSO) )INDSX( I , J ) =1 
I F ( (OUT .GT.PADSQ ) . AND . ( IN.GT.RADSQ ))IND5X( I , J > =0 
I F ( ( OUT. GT. RADSO ) . AND. ( I N. LE . RADSO ) )IND2X( I, J ) = 2 

41 CONTINUE 

42 CONTINUE 



QUADRANT II CHECK 
DO 44 J=1 , 21 
DO 43 1=1. 21 

OUT=(V*TL-( 1-1 )*DELT AX )**2+ ( V*TL- ( J-l ) *DELTAX )**2 
IN = ( V*TL-I*DELTAX)**2+( V*TL-J*DELTAX)**2 
RADSO= (V*TL )**2 

I F ( (OUT. LE.RADSO). AND. (IN. LE.RADSO) ) INDEX ( I,J) = 1 
I F < ( OUT. GT.PADSQ). AND. (IN.GT.RADSQ) ) I NDEX( I ,J)=0 
I F ( ( OUT . GT . RADSO ). AND. ( IN. LE .RADS Q ) )INDEX( I,J) = 2 

43 CONTINUE 

44 CONTINUE 



QUADRANT III 
DO 46 J=l, 21 
DO 45 1=21,41 

OUT=( V*TL-( J-l)*DELTAX)**2+( I*DELTAX-V*TL) **2 
IN=(V*TL-J7DELTAX )**2+( ( 1-1 )*DELTAX-V*TL )**2 
RADSQ=( V*TL)-?*2 

I F ( (OUT. LE.RADSO ). AND. ( IN. LE.RADSO) ) I ND5 X ( I , J ) =1 
I F ( ( OUT. GT. RADS 0). AND. ( IN.GT.RADSQ) ) I ND6X ( 1 , J ) = 0 
IF ( (OUT. GT. RAD SO). AND. (1 N. LE.RADSO) )INDEX( I , J ) =2 

45 CONTINUE 

46 CONTINUE 



QUADRANT IV CHECK 
DO 48 J=21 ,41 
DO 47 1=21,41 

OUT=( V*TL-I» DE LT AX )^*2 + ( V*TL -J*DSLT AX )** 2 

IN =( V*TL-( 1-1 J^DELTAXJf *2+( V*TL-( J-l )*DELTAX )**2 

RADSQ= (V-TL )**2 

IF( (OUT. LE.RADSO). AND. (IN. LE.RADSO) ) INDEX ( I,J)=1 
I F ( ( OUT. GT. RAD SO). AND. (IN.GT.RADSQ) ) I NDS X ( I , J ) =0 
I F ( ( OUT. GT. RADS Q). AND. ( IN . L E .RADSO ) ) INDEX( I, J ) = 2 

47 CONTINUE 

48 CONTINUE 



NOW TO SIMULATE THE PARTICLE MOTION UNDER THE 
SPECIFIED RULES. THIS IS THE MAIN DO LOOP IN THE 
PROGRAM. FIRST A TIME UNTIL THE NEXT COURSE CHANGE 
IS GENERATED. THEN A COURSE TO BE STEERRED IS GEN- 
ERATED. XX AND YY ARE THE TWO COMPONENTS OF POSITION 
ADDED DUE TO A COURSE AND TIM? SEGMENT. X AND Y ARE 
THE UPDATED POSITION c ROM DATUM. AT TIME-LATE THE 
UPDATED POSITION IS CONVERTED TO A CELL POSITION AND 
THAT CELL VALUE IS INCREMENTED BY ONE. 
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DO 3 J = l, 16810 

X=0.0 

Y=0.0 

CLOC K=0 • 0 

C=0 .0 

TR=TL-C LOCK 

CALL RANDU( I X , I Y ,R ) 

IX = IY 

CALL CORSE ( R » C ) 

THETA=C 
GO TO 66 
1 TR=TL-C LOCK 

CALL RANDU(IX,IY,R) 

I X = I Y 

CALL TPMCRS (R »C ) 

66 CALL RANDU< I X,I Y,R) 

IX = IY 

CALL T RKTM (FT » R, T ) 

TT=AMI Ml ( T , TR ) 

CALL MOV E ( V , TT ,C»XX,YY) 

CLOCK=CLOCK+T 

X=X+XX 

Y=Y+YY 

IF(CLOCK.LT.TL)GO TO 1 

XPRI Mc*= X^COS ( THETA ) +Y*SIN( THETA) 

YPRI ME=Y*COS( THETA ) -X*S IN ( THETA ) 

X=XPRI ME 
Y=YPR I M 5 

CALL RANDU ( IX, IY,R) 

I X=I Y 
C0RR=1 .0 

CALL HISTO(K,DELTAX,R,CORR,X, I ) 

L = I 

CALL RAMDU ( I X, I Y ,R ) 

IX = I Y 
CORR=- 1.0 

CALL HISTO(K,DELTAX,R,CORR,Y,I ) 

M = I 

C THE PARTICULAR CELL I,J IS DETERMINED BY TWO SUCCESS- 

C IVE CALLS OF HISTO. IT IS POSSIBLE A POSITION X,Y 

C COULD BE DETERMINED TO Be IN AN INFEASIBLE CELL. 

C IF THIS HAPPENS THEN RESOLV WILL BE CALLED TO CHANGE 

C THAT CELL ASSIGNMENT IN A PRESCRIBED MANNER TO A 

C FEASIBLE CELL. 

IF(IND5X(M,L ).5C.0)G0 TO 31 
30 PP ( M , L ) =PP ( M, L )+l 
GO TO 3 

21 CALL RE SOL V( L » M , I , N ) 

NC0UNT=NC0UNT+1 

L = I 

M=N 

GO TO 30 
3 CONTINUE 
C 

NUPC5L =0 
NUPEDC=0 
CRITS2=0.0 
DO 52 1=1 ,41 
DO 51 J = 1 » 41 

I F ( P P ( I , J) .EO.OIGO TO 51 

NUPCSL=NUPCcL+l 

I F( I NDE X ( I , J ) • NE • 2) GO TO 51 

NUPEDC=NUP2DC+1 

CRI TE2 = C RI T-2 + PP (I , J) /1 6 81 0.0 

51 CONTINUE 

52 CONTINUE 

CR I TE3 = NUP fr DC 
RNUPED=M UPEDC 
RNU P C 2 = NU P r E L 
CRITEl=RNUP£D/RNUPC£ 
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NUPINC=NUPCEL-NUPEDC 
212 WR!T=(6,203) 

203 FORMAT ( * 1 * ,60X , ' PP PAGE 1') 

DO 205 1=1.41 

WRITE (6 . 204 ) ( PP( I , J ) ,J=1, 15 ) 

204 FORMAT ( ' • ,1517) 

205 CONTINUE 
WRITE (6,206) 

206 FORMAT ( • 1 • ,60X , • PP PAGE 2') 

DO 208 1=1,41 

WRITE(6,207) (PP(I,J),J=15,29) 

207 FORMAT ( • • , 151 7) 

208 CONTINUE 
WRIT? (6,209) 

209 FORMAT ( • 1 • » 60X . • PP PAGE 3') 

DO 211 1=1 ,41 

WRI TE(6,210) (PP( I ,J) ,J=29,41) 

210 FORMAT ( • • , 1317) 

211 CONTINUE 

NOW THAT THE HISTOGRAM IS BUILT ITS ELEMENTS ARE 
ORDERED SO THE HIGHEST VALUES ARE IN THE LOWEST NUM- 
BERED CELLS. 



1000 IC=0 

DO 1010 1=1,1680 

IF(XP( I ) .G5.XP ( I+l ) ) GO TO 1010 

MAX=XP( I+l > 

XP ( I +1 )=XP ( I ) 

XP ( I ) =MAX 
IC = 10 

1010 CONTINUE 

IF(IC.NE.O) GO TO 1000 
DO 2C00 1=1,41 
DO 2010 J=1 , 41 
OQ ( I , J) =PP( J,I ) 

2010 CONTINUE 
2000 CONTINUE 

DO 1030 1=1,41 
DO 1020 J= 1 , 41 . 

PP(I , J)=00(I,J ) 

1020 CGNTINUE 
1030 CONTINUE 



AFTER ORDERING THE CELL VALUES THEY ARE ACCIJMLATED SO 
EACH CELL CONTAINS THE SUM OF ITS OWN VALUE AND ALL 
THE VALUES PRECEDING IT. 



CUMSUM=0 
DO 11 1=1,41 
DO 10 J = 1 , 41 
CUMSUM=CUMSUM+PP(I , J) 
PP ( I , J)=CUMSUM 

10 CONTINUE 

11 CONTINUE 



NUMSXC=0 
NUMI NC = 0 
NUMEDC=0 
DO 33 1=1 ,41 
DO 32 J = 1 , 4 1 

IF ( I N05X ( I , J ) . SQ. 0 ) NUMEXC=N U W EXC + 1 
IF ( I NO E X ( I , J ) .50.1 ) NU M I NC = NU M I N C + 1 
IF ( I NDE X( I,J).EO. 2)NUMEDC=NUMEDC + 1 

32 CONTINUE 

33 CONTINUE 

numcel=numinc+numedc 

totcel=numexc+numedc+numinc 
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NOW TO LOAO THE FIRST 900 CELL VALUES INTO THE PLOT- 
TING VECTOR EY. THE CUMULATIVE NUMBER OF LOOKS OR 
KILL POWER WILL BE PLOTTED ON THE X AXIS. 



EEX( 1 ) =0.0 
EEX ( 2 ) = 1 .0 
EEY( 1) =0.0 
EEY( 2)=1.0 



DO 15 1=1,16 
DO 14 J=l,41 
N=(I -1)*41 + J 
EY! N )=PP ( I » J J/16810.0 
EX (N ) =N/13 85.0 

14 CONTINUE 

15 CONTINUE 
DO 16 J=1 ,36 
1 = 17 

N= ( I -1 )*41+J 
EY ( N ) =P P ( I ,J )/16810.0 
EX(N)=N/1385.0 

16 CONTINUE 
C 

c 

DO 17 J=37,41 
1=17 

N=( I-l)*41+J 
L=N-69 2 

EYY ( L ) =PP( I , JJ/16810.0 
EXX(L)=N/1385.0 

17 CONTINUE 

DO 19 1=18,33 
DO 18 J= 1, A 1 
N=U -1 )*41+J 
L=N-692 

EYY (L )=PP ( I » J ) /1 68 10.0 
EXX(L)=N/1 385.0 

18 CONTINUE 

19 CONTINUE 

DO 20 J=1 ,31 
1=34 

N=( I-l)*41+J 
L=N-692 

EYY ( l. ) = PP( I ,J) /16810.0 
EXX(L)=N/1385.0 

20 CONTINUE 

WRITF( 6, 606) NCOUNT 

606 FORMAT!//, 25X, 'NCOUNT EQUALS ',116) 

WRITE (6,100) NUM2 XC , NUM5DC , NUMI NC , TOTCSL 

100 FORMAT ( / * 10X, 'NUMEXC= ' , II 6, 5X, 'NUMEDC = ' , II 6, 5X, ' NUMIN 
XC = ' ,1I6,5X,'T0TCEL=',1I6) 

WRITE! 6, 101 ) CRI T2 3 ,CRI TE1 »CRI TE2 

101 FORMAT (//, 10X, • NUMBER OF POSITIVE EDGE CELLS IS ',11 
X6 , / / , 10 X , ' THE PERCENTAGE OF POSITIVE CELLS WHICH ARE E 
XDGE CELLS IS • , IF A. 3, // , 10X , • THE AMOUNT OF PRQ8A8ILIT 
XY IN THE EDGE CELLS IS ',1F4.3) 

CALL DRAW ( 2 ,EE X ,EE Y ,1 ,0 , LABEL , I T I TL5 ,0 ,0 , 0 ,0 ,0 , 0 , 9, 9, 
XO* L AST ) 

CALL DRAW (692 , EX, EY , 2 ,0 , LABEL , IT ITL 5 , 0 , 0, 0, 0, 0 , 0 , 9 , 9, 
XO , LA ST ) 

CALL DRAW (692,EXX,EYY,3,0,LABEL,ITITLE,0,0,0,0,0,0,9, 
X9 , 0 , LAST ) 

STOP 

END 
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SUBROUTINE RANDU ( I X , I Y , YFL ) 

IY=IX*65539 

IF(IY) 5,6,6 

5 IY=IY+21 47 483647+1 

6 YFL= IY 

YFL=YFL*.4656613E-9 

RETURN 

END 



SUBROUTINE TRKTM (E T ,R , T) 
RA=1.C-R 
Tl=-E T* ALOG ( R) 
T2=-5T*’AL0G(RA ) 

T = T1 

RETURN 

END 



SUBROIJTI N5 MOVE( V2L,TI NE , ANG ,X ,Y ) 
X=VEL*TIME*COS(ANG) 

Y=VEL' i: TI ME^SIN ( ANG ) 

RETURN 

END 



SUBROUTINE HI S TO ( K , D ELTAX , R , COP.R , X , I ) 

UP3ND= ( ( 2*K+l)*0EUTAX)/2.0 

L=2*K+1 

XX=( X+(CORR*UP8ND) ) /DELTAX 

XX=CORRT ; XX 

XI =AMOD ( XX ,1 .0 ) 

M=INT( XX) 

I F ( M . EQ .0 ) GO TO 1 
I F ( M.EQ.UGO TO 2 
IF( Xl.Ns.O.OJGO TO 3 
I F ( R ♦ GT • 0 • 5 ) GO TO 3 
I =M 

GO TO 4 
1 1 = 1 
GO TO 4 
2 1 = 1 
GO TO 4 

3 I =M+ 1 

4 RETURN 
END 



SUBROUTINE RE SOU V< I , J , L , M) 
IFd.GT .21) GO TO 2 
IF(J.GT. 21) GO TO 1 
U= I + 1 
R = J + l 
GO TO 4 

1 L=I+1 
M = J-l 
GO TO 4 

2 I F ( J .GT.21JG0 TO 3 
U = I-1 

R = J+ 1 



R = J-1 
4 RETURN 
END 
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SUBROUTINE CORSE(R,C) 

PI=3. 1415926 

A=0. 0 

B=2.0*PI 

F=2.0*PI 

DELTAC=R*< B-A ) +A 
C=C+DcLT AC 
C=AMOD ( C » F ) 

RETURN 

END 



SUBROUTINE TC0RS5(R,C) 

R IS THE UNIFORM VARIATE USED TO GENERATE THE NEXT 
COURSE CHANGE. 

C IS THE PRESENT COURSE AND WILL BE RETURNED AFTER 
ADDING A NEW HEADING CHANGE. 

DELTAC IS THE AMOUNT OF THE NEW HEADING CHANGE 
S=0. 20 

S IS THE TEEPEE SLOPE AND IS BETWEEN 0 AMD 
2/(PI**2) = 0.202. NOTE , S CANNOT EQUAL 0. 

PI=3. 1415927 

IFIO.O.LE.R.AND.R. LE. 0.25) GO TO 1 
IF(0.25.LT.R.AND.R.LE.0.5)G0 TO 2 
IF(0.5.LT.R.AND.R.LE.0.75)G0 TO 3 
IF(0.75.LT.R.AND.R. LE. 1 • 0) GO TO 4 



A=S/2. 0 

B=(1.0/(2.0*PI ) )-< S*P 1/4.0) 

D=-l . O rt R 
E = B^ ; / : 2-4. 0*A*D 

D£LTAC= (-1.0*B+SQRT(E ) )/( 2.0*A ) 
GO TO 5 



3 R PRI ME = R— 0 .5 
A= S/2. 0 

B= < l.O/I 2.0*PI ) )-< S^P 1/4.0) 

D=-l ,0*R PRIME 
£=8**2- 4. 0*A* : D 

DELTAC= <-1.0*B+SQRT(E ) )/( 2.0*A) 

DELTAC=DELTAC+ D 1 

GO TO 5 



C 

C 



2 A= ( -1. O yt S) /2.0 

B=( 1.0 / ( 2.0* PI ) ) +( 3.0*S*PI /( 4.0) ) 
D=-l .0=M (S*PI**2 )/4 .0 ) -R 
E= 2- 4 . 0*A *D 

DELTAC=(-l.C 3 ‘B+SQRT(E))/( 2.0*A) 

GO TO 5 



4 RPRI ME=R-Q .5 
A = ( - 1 . 0* S ) /2.0 

B=(1.0/ (2.0*P I ) ) + ( 3.0*S*PI/< 4.0) ) 
D=-l .0*( <S*P!**2 )/4.0 ) -RPRI ME 
E=B**2-4.0*A*D 
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DELT AC= ( -1 .0*3 +S QRT (E ) )/(2.0*A) 

DELTAC=DELTAC+PI 

GO TO 5 



C 

C 

C 

C 



5 C=C+DELT AC 
F=2.C*PI 
C=AMOD( C t F ) 



RETURN 

END 



SUBROUTINE PCORSE(R.C) 

PI =3. 1415926 
F=2.0*PI 
S=1.0/ ( °I**2) 

IF(R.GT.O. 5 ) GO TO 1 

A=2.0*P/S 

DELT AC = S QRT ( A) 

GO TO 3 

1 A=- 1 ♦ 0 s " S/2 . 0 

B=(1.0/(2.0*PI ) )+(3.0*PI*S/2.0) 

D=-l . 0* S*(P!**2)-R 
E=(8**2)-4.0*A*D 
D5LTC1=(-B+SQRT(=) )/(2.0*A) 

D£ LTC 2= ( -8 -SORT ( E ) ) /(2.0*A) 
IP(PI.LT.DELTC1.AMD.DELTC1.LE.F)G0 TO 2 
CELT AC=DELTC2 
GO TO 3 

2 DELTAC=DELTC1 

3 C=C+DELTAC 
C=AMOD(C ,F ) 

RETURN 

END 



SUBROUTINE LRCQRS(R.C) 
PI =3. 1415926 
F = 2.0* PI 

IF(R.GT.0.5)G0 TO 1 
C=C+ ( PI/2 .0 ) 

GO TO 2 

1 C=C + ( <3.0*PI )/2.0) 

2 C=AMQD(C,F) 

RETURN 

END 



SUBROUTINE TRNCRS(R.C) 
PI=3. 1415926 
F=2.0*PI 
A=PT / 6. 0 
B=(11.0*PI ) / 6 . 0 
DELTAC=R*( B-A) +A 
C=C+DEL T AC 
C= AMOD ( C » F ) 

RETURN 

END 
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